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Abstract: The dynamical systems methods are used to study evolution of the 
polymerised scalar field cosmologies with the cosmological constant. We have found 
all evolutional paths admissible for all initial conditions on the two-dimensional phase 
space. We have shown that the cyclic solutions are generic. The exact solution 
for polymerised cosmology is also obtained. Two basic cases are investigated, the 
polymerised scalar field and the polymerised gravitational and scalar field part. In 
the former the division on the cyclic and non-cyclic behaviour is established following 
the sign of the cosmological constant. The value of the cosmological constant is upper 
bounded purely from the dynamical setting. 
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In the last years the procedure of background independent quantisation has 
been applied to the mini-superspace cosmological models. This approach is known 
as Loop Quantum Cosmology ]IJ. In this theory the gravitational part is quantised 
in the background independent way (polymer quantisation) while the matter part 
remains classical. The significant result of this theory is resolution of the cosmological 
singularity problem. This was strictly proved only for the models with the free scalar 
field. However the considerations based on the "effective" equation show that it is 
the case also for some models with a potential function ||. 

Recently the procedure of polymer quantisation has been applied also to the mat- 
ter part. In their paper Hossain et al. || constructed the model with polymerised free 
scalar field. In this model the gravitational part is treated classically. They showed 
that the evolution is nonsingular because the energy density of the polymerised is 
always finite. Moreover the inflationary phase emerges from this model. 

The natural continuation of these studies would be construction of the fully quan- 
tum polymerised cosmological model. In such a model both the gravity and matter 
part would be quantised in the background independent way. However Hossain et al. 
[[J indicated the possible problems with construction of the proper algebra for poly- 
mer variables. It is due to the fact that matter part variables involve gravitational 
degrees of freedom. 

In this paper we construct the classical model with the known effects of polymeri- 
sation. The model is classical therefore the problems indicated on the quantum level 
do not occur. The effects of polymerisation are introduced by the phenomenological 
factors in the classical Hamiltonian. We do not state that the model reproduces the 
results which could be obtained in the fully quantum approach. Perhaps in some 
regimes, where quantum correlations between matter and gravity are strong, this 
approximation does not work. In this approach we assume that the gravity factors 
in the definitions of polymer variables of the scalar field are the mean values. In 
particular in the definition of the polymerised momentum operator we have 



where p is a geometric area operator. Based on this assumption we derive "effective" 
equations of motion for the system. We set the FRW symmetry of the gravitational 
part and assume matter content to be a free scalar field together with cosmological 
constant. The classical phase space of this model is parametrised by the canonical 
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variables (c,p, <fi, tv^) which fulfil the Poisson brackets {c,p} 
The classical Hamiltonian for this model is given by 



f^and {0,tt4 = 1. 
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The effect of polymerisation can be now introduced by the following replacement in 
the above Hamiltonian 

p 3/2 f A 



— sin^-^J. (4) 

Based on this we obtain the phenomenological Hamiltonian 7i p hen which depends on 
the unknown parameters \x and A. The dimension of \p] = L and [A] = L 2 therefore 
they can be interpreted respectively as a length and area scale of polymerisation. 
We do not assume here any relation between \x and A. The polymerisation of gravity 
and matter part can have quite different energy scales. In the rest of the paper we 
will investigate two models, first one, based on the polymerised quantisation of the 
free scalar field, and second, where the effects of polymerisations are included both 
in the gravitational sector and the free scalar field. 

We concentrate on the topological structure of the phase space and different 
evolutional paths of the investigated systems which can be put in the simple form 
of a dynamical system of the Newtonian type |Q. Within a large class of solutions 
the special role is played by the cyclic solutions f| || which are generic (although 
structurally unstable). They are located around a centre type critical points. 

For the first model we investigate the system described by the following Hamil- 
tonian 

H phen = ~^VP+ ^2 Sin 2 (-i^y/ 2 + V 2 (5) 

where k = 8nG. Based on the Hamilton equation = {/, H p hen} together with the 
Hamiltonian constraint Ti. p \ lcn ~ we can derive the modified Friedman equation 

where the energy density is 

Here tt^ = const, and is a parameter of the theory. Moreover the equation for the 
evolution of the field <fi is given by 



# 1 . / A \ 
d^ = 2A Sm 



In order to perform dynamical analysis of the investigated system first we make the 
following change of the dynamical variable 

A 

x ~ ^v^ 71 ^' 
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next the time reparameterization 



;2A 2 1 , 
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then the modified Friedman equation can be put in the following form 

H 1 = -x' 2 --sm 2 (x) = -\ 2 . (9) 
2 2 k K ' 

The resulting dynamical system is of the Newtonian type 



y 



av{x) (10) 



where the potential function is 



V(x) = -^sm 2 (x), (11) 



and the term proportional to the cosmological constant A plays the role of constant 
energy level. 

At the finite domain of the phase space any system of the Newtonian type has 
the critical points of a saddle or a centre type only. The character of a critical point 
is determined from the characteristic equation of the linearization matrix, namely 

f= d 2 V(x) 



dx 2 

x 

calculated at the critical point x Q . If I 2 > we have a saddle type critical point, in 
opposite case if I 2 < we have a centre type critical point. 

For investigated system (|l(J) coordinates of the critical points are (xo,yo) = 
(§7r,0) , k = 0,1,2,..., and their character can be identified as: a saddle type 
critical points x$ = kn and a centre type critical points xq = (2k + l)n where 
k = 0,1,2, .... 

The phase portrait for this system we present in Figure [I]. The dashed trajecto- 
ries correspond to the vanishing cosmological constant case. Clearly there is a family 
of periodic solutions which are only possible for the negative cosmological constant. 

Equation @ can be easily integrated for different special cases. 

• A = 0: 

r + const. = In (tan (12) 



A > 0: 




r + const. = ; : F [ x\ - oA ) (13) 
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where F(-\-) is an incomplete elliptic integral of the first kind 



- 3 - 



• for A < we have series of periodic solutions where the period is given by 
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Figure 1: The phase portrait for the first model under consideration where the scalar 
field is polymerised. All trajectories lies in the physical region. The dotted lines connecting 
saddle type critical points denote trajectories for the vanishing cosmological constant A = 0. 
Regions around centre type critical points are occupied by the trajectories for negative 
values of the cosmological constant A < 0. The rest of the phase space in occupied by the 
trajectories for A > 0. 
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Figure 2: The evolution of the scale factor (black) and the Hubble function (red) for a 
negative value of cosmological constant in the first model under considerations. This is a 
sample of a cyclic evolution in the model. The bounce and the recollapse take place for 
vanishing total energy density. This is a specific feature of the model where the free scalar 
field is polymerised. 
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For the second model we investigate effective dynamics described by the Hamil- 
tonian constraint where the gravitational and free scalar field parts are both poly- 
merised 

- (^> 3/2 + v <«» 

The modified Friedmann equation is in the form 

* - h - j) < 16 > 

where the total energy density 

1 • 2 M A A 

P T = 0T2 S111 (375^) +- ( 17 ) 



2A 2 Vp 3/2 
and the critical density 

p c = -4-n- 

Note that from (|16D and (|P7D follows the upper limit on the positive cosmological 
constant 

A < ^Pcrit- 

Similar dynamical restriction on the value of the cosmological constant was also 
found in the flat model with a scalar field and holonomy corrections of LQG |7|, ||. 
The modified Friedmann equation can be put in the following form (dividing by 

Pt) 

1 -- 2 ^- = - 1 -s^(x)+ 1 -a-h> (19) 
2 x z Skpt 2 2k 

where a = 2X 2 p c . After the following time reparameterization 



dt=J 2 ^ 1 -da (20) 
y 3k pt x 

we receive 

H 2 = \x' 2 -\(a- sin 2 (xj) = --A 2 (21) 

In Figure [3] we present the phase portraits for different values of the parameter a. 

Note that the parameter a depends on both scales of polymerisation // and A. If 
we for now assume that for the gravitational part we have standard Jx scheme then 
p 2 = A = 2V3irjl 2 P and 

6 A 2 V3 x2 
k^ 2 p z 87r 2 7 d /p 

From the phase space analysis Fig. [| we can notice that for vanishing cosmological 
constant the qualitative behaviour is different for a > 1 and a < 1. Therefore the 
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a = 1 is a boundary value for which we can calculate from the previous expression 
the boundary value of A parameter, namely 



'5* <*> 



Next we can present investigated system (21) in the form of a Newtonian type 



dynamical system with the potential function given by 

V(x) = --{a- sin 2 (x)) . (24) 

As in the former case the organisation of the phase space depends only on the form 
of the potential function. The location of the critical points is the same as in the 
former case Xq = |7r, k = 0, 1,2, . . . , but their character is reversed, namely, for 
xo = k% we have a centre type critical points and for xq = (2k + 1)tt we have a 
saddle type critical points. On the Fig. |3| we present the phase space portraits for 
the system for different values of the a parameter. 
Now Eq. (|2l"l) can be easy integrated 



dx\ 2 2 f A 
da / \° n 



X 2 - sin 2 (x) (25) 



o + const. = - =F x\— x — (26) 

^2 ( Pe - A) \2 \ 2(p c --)P 

where F(-\-) is an incomplete elliptic integral of the first kind. 
The period of the periodic solution is given by 

T a = ; 4 F I arcsin I W2 ( p c - — 1 A 2 I \ -, | (27) 



V 2 (P* - I) A2 




The unexpected feature which distinguishes this model from those meet in the 
standard loop quantum cosmology is that for periodic solutions the bounce and rec- 
ollapse take place for the total energy density equal the critical density. This happens 
due to presence of the periodic function sin 2 (x) in the total energy density. On the 
Fig. |] we present the cosmological time evolution of three quantities proportional 
to the Hubble function H 2 , a difference between the total energy density and the 
critical density p c — pt and finally the scale factor a. Initial conditions was chosen 
in such a way that at t = the universe was in its maximal expansion phase. 

The main conclusion of this letter is that the cyclic behaviour is generic feature 
of polymerised cosmology. Two version of polymerisation has been considered in 
the model with cosmological constant. First, the polymerised scalar field only was 
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studied, and second, both gravitation part and scalar field have been polymerised. 
The cyclic trajectories are represented in the phase space by closed phase curves 
around a centre type critical point. We have found the expression for their periods 
as well as exact solutions for trajectories. It is interesting that the cycle is asymmetric 
in original dynamical variables. Because of the reparameterization this information 
is hidden. 

The dynamics of the model also offers the upper limitation on the cosmological 
constant value. 
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Figure 3: The phase diagrams for the second model under consideration, with both 
gravitational and scalar field parts polymerised. From top to bottom: a = h, a = 1, a = ^. 
On the line x = p is infinite we have a singularity. The dashed lines denotes trajectories for 
the vanishing cosmological constant. In the regions between them trajectories for positive 
cosmological constant are located. The shaded regions are unphysical. Every trajectory 
initially located in the physical domain which tend to the unphysical region reach H 2 = 
after an infinite period of the cosmological time. Note that in all cases the trajectories for 
the zero cosmological constant are tangent to the border of the unphysical region for the 
same value of the coordinate x. The centres have coordinates x = kir, k = 1,2,... and 
saddle have coordinates x = (2k + l)f , k = 0, 1, — 






- 9 - 



